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Abstract 

The class of deterministic 'Daphnia' models treated by Diekmann et al. (J Math Biol 61: 277-318, 
2010) has a long history going back to Nisbet and Gurney (Theor Pop Biol 23: 114-135, 1983) and 
Diekmann et al. (Nieuw Archief voor Wiskunde 4: 82-109, 1984). In this note, we formulate the 
individual based models (IBM) supposedly underlying those deterministic models. The models treat 
the interaction between a general size-structured consumer population ('Daphnia') and an unstruc- 
tured resource ('algae'). The discrete, size and age-structured Daphnia population changes through 
births and deaths of its individuals and throught their aging and growth. The birth and death rates 
depend on the sizes of the individuals and on the concentration of the algae. The latter is supposed 
to be a continuous variable with a deterministic dynamics that depends on the Daphnia population. 
In this model setting we prove that when the Daphnia population is large, the stochastic differential 
equation describing the IBM can be approximated by the delay equation featured in (Diekmann et al., 
I.e.). 

Keywords: Birth and death process; age and size-structured populations; stochastic interacting particle 
systems; piecewise deterministic motion; large population limits. 
AMS codes: 92D40; 60J80; 60K35; 60F99. 

1 Introduction 

The theory of physiologically structured populations as developed in works by Diekmann and co-authors 
[2~il [7J [5J HD] derives its motivation from its supposed ability to link population level phenomena to 
specific mechanisms in and around individual organisms. Yet, those individuals do not figure as such in 
the models, which treat the spatial concentrations of those individuals as continua. The supposition has 
always been that the proposed deterministic frameworks would in principle be derivable from individual- 
based stochastic models (c.f. [2S]), but full derivations of this sort so far have only been done for the 
special cases of finite i-state models (i from individual; leading to ODEs in the large number limit, 
e.g.[ini[H]) and age-based models (e.g. [23 [Ml 1301 EES El EES]). Moreover, it is possible to reinterpret 
the results in [T] as pertaining to general structured models with only jump transitions. Although some 
age-structured models are phrased in terms of i-state variables that change with time in a fixed manner, 
uninfluenced by the environment, these can truly represent but a small fraction of the rich variety of 
physiological mechanisms seen in nature (e.g. [41 1241 chapters I and III]). For example, for the majority 
of species, size with a growth rate that depends on the environment is a far more important determinant of 
an individual's population dynamical behaviour than is age. Except in the physiologically well buffered 
homeotherms (on which we anthropocentrically are inclined to focus) reproduction tends to be under 
a strong influence of past food availabilities, with reproductive size being reached far earlier when an 
individual has encountered good than when it has encountered bad feeding conditions. In, for example, 
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arthropods or fishes, maturing to the reproductive states often largely depends on how much food they 
have encountered but little on their age. However, the toolbox for proving the appropriate law of large 
number results that was developed in [Ml [29j [30l [13] needs to be extended a bit to deal with size-based 
models, referred to as 'Daphnia' models by [10], let alone to the even more general models considered in 
[2"3l [231 [SJ El IE] - In this paper we, as a birthday present to Odo Diekmann, intend to remedy this lack, 
although so far only for the 'Daphnia' models that appear as central example in his work (c.f. |10|). 

Diekmann et al. |10j lay out a general framework for investigating the interaction between a general 
age or size-structured consumer population (referred to as 'Daphnia') and an unstructured resource (re- 
ferred to as 'algae'), a class of problems special forms of which were first considered in [26l [15] and [51128]. 
Diekmann et al. (I.e.) show that stability properties and bifurcation phenomena can be understood in 
terms of solutions of a system of two delay equations that are analysed using results in [HIGH- In this note, 
we derive their model from a microscopic description: starting from a stochastic, age- and size-structured, 
individual based model (IBM) for the Daphnia population we recover the equations of [TU] . 
In our study, the Daphnia population is discrete and stochastic while the algal population is continuous. 
We treat the Daphnia population as a point measure on a space spanned by size and age and from 
this platform generalize the microscopic construction given in [30) to arrive at a stochastic differential 
equation (SDE in the sequel) driven by Poisson point processes that gives a pathwise description of the 
population evolution. The main novelty is that the individual growth rate also depends on the population 
as a whole, instead of only on an individuals' characteristics, since the algal concentration and thus each 
individuals' resource access is affected by the entire Daphnia population. 

We then, in the wake of [371 I2H1 1301 1131 [HI HE] > provide a law of large numbers that allows approximating 
the dynamics of the individual-based size- and age-structured model by means of a partial differential 
equation (PDE) when the volume containing the Daphnia and algae as well as the population sizes are 
large. In |23j similar limit theorems are used as a basis for obtaining adaptive dynamics approximations 
for the evolution of hereditary traits in age-structured populations, while [11] gives informed conjectures 
about such approximations for the more encompassing model classes treated in [5]. 

Finally we study the limiting equations to recover the macroscopic (deterministic) system considered in 



Notation: We will denote the set of finite measures on endowed with the weak convergence topology 
as Mp(R^_). For fi £ _M_f(R^) and a real measurable bounded function / on this space, we set (/x, /) = 
Jr2 f(£,a)/j,(d£,da). The set of bounded real functions of class C 1 on E with bounded derivatives is 

denoted as C^(E,R). 

If X is a process indexed by time, then we will denote the value of X at time t as either X(t) or X t (the 
latter to avoid formulas becoming cluttered with too many brackets). 

2 Individual-based Daphnia model 

2.1 Model specification 

Our Daphnia population consists of discrete individuals living in continuous time, differentiated by an 
age a £ R+ and a size £ € R+. Individuals are given labels i g N* = {1,2,., .}, with the individuals 
present in the population at t = bearing numbers in an order opposite to that of their ages, followed 
by the individuals born after t = in the order of their appearance in the population. We denote as 
It G N* the total number of individuals that were present at time or were born between time and 
time t. Vt C {1, . . . , It} denotes the set of individual alive at time t. We can then represent the Daphnia 
population as a point measure on the state space N* x R + x R + : 



where Q and a\ are the size and age of individual i at time t. We denote as Z t (d£, da) = Z t (N* x d£ x da) 
the marginal measure of Z t on R+ ; Vt equals the support of the marginal measure of Z t on N* . Size and 
age are related as follows. 

An individual's age is equal to t — to where t is the current time and to its birth time. Individuals grow 
up in an environment that is characterized at time t by the algal concentration S(t). All individuals are 
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born with the same size £o (to keep things simple) and an individual with size £ at time t grows at speed 
<?(£, S{t)), so that the size of an individual aged a born at time t$ is: 

£Mo)=£o+ / g(Z(a;t ),S(t +a))da. (2) 
Jo 

The growth rates depend on the other individuals in the population. The latter point is the novelty of the 
IBM presented here. (We could have taken hereditary traits as well as further age-like i-state variables 
on board as in [30j . but have decided not do this in order not to unduly complicate the story). 
Reproduction is asexual. The birth and death rates of an individual with size £ and age a at time t are 
/?(£, a, S(t j) respectively /x(£, a, S(t)). 

An individual with size £ depletes the food density at rate 3^7 (£, S(t)) and the food density replenishes 
with rate f(S(t)) so that the food concentration S(t) evolves according to 

f (*) =f(St) - ± [ 7(Z,St)Z t {de,da) = f(S t ) - i rt&St). (3) 

For the biological justification think of the Daphnia population as living in a container of size K, so that 
the Daphnia density is with \Z t \ :— Card(Vt) = (Z t , l) the number of individuals. 

Assumption 2.1. In the sequel, we assume that the growth speed g(£,S), the rates f(S) and j(£,S) 
are continuous bounded functions and that 

(i) the birth rate /3(£, a, S) is piecewise continuous and bounded by /3. 

(ii) the death rate a, S) is continuous and there exists a function /i(a) and a constant A G (0, +oo] 
such that V(£, a, 5) G K^_, a, 5) > /i(a) and Jjj 4 fi(a)da — +oo. 

(iii) 5 is Lipschitz continuous with respect to £, uniformly in S on compact intervals of R+, and bounded 

by g- 

(iv) / and 7 are uniformly Lipschitz continuous with respect to S uniformly in £ on compact intervals of 
R+. 

The assumption (i) on the birth rate ensures that in a short time interval a single individual can not 
beget too many young: intervals between births are stochastically lower bounded by exponential random 
variables with rate /3. The assumption (ii) on the death rate implies that individuals a.s. have lifetimes 
bounded by A. Finally the assumptions (iii) and (iv) ensure that there exist unique continuous solutions 
to (0) and §3§ as long as the number of individuals \Z t \ remains finite, the latter being guaranteed by the 
fact that \Z t \ is stochastically bounded by a pure birth process with birth rate /3. 

For $ G Cb(M 2 ,IR) and ip G Cl(R%,R), we denote by <J> the function on M f (l 2 + ) x M + defined by 
$ V (Z,S) = $>((Z,ip),S). From the description of the population dynamics, it follows that the process 
(Z(t), S{t))tem. + is characterized by the infinitesimal generator L operating on the functions $ ¥ ,: 

L%(Z, S) - d^Z, S) (Z, g(., S),dtf(., .) + d a f(,, .)) 

+(Z, M (., ., S) tp) - <p(., .), S) - %{Z, S)) ) 

+(Z, /?(.,., S) (#«Z, ip) + ^o,0), S) - %(Z, 5)) ) 
+d 2 %(Z,S)(f(S)-(Z,± 1 (.,S))) (4) 

The first term describes the aging and growth of the living individuals of the population. The second 
and third terms represent the demography of the population (deaths and births) . The fourth term cor- 
responds to the variation of the food. 

In the next Subsection we introduce the pathwise construction of an IBM with the described dynamics 
and give an SDE driven by a Poisson point process that admits (j4]) as infinitesimal generator. This is 
useful for simulations and for deriving moment conditions and large population approximations (e.g. [T]). 
The evolution is piecewise deterministic: The size of the population is modified at birth or death events. 
Between these, conditionally on the structure of the population after the last event, the growth of the 
individuals and the food dynamics are deterministic. 
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2.2 Construction of the IBM and a useful SDE 



Let us start with some heuristics. Consider at time t a population given by Z t and food concentration 
S(t). If no birth or death event occurs between time t and t + s, then Vt+ S — Vt- The sizes £ l (t + s) for 
i £ Vt at time t + s and the food concentration S(t + s) are obtained by solving: 

C(t + s)=e(t) + 9(e(r),S(r))dr (5) 

S(t + s) = S(t) + f +S (f(S(r)) - 1 £ 7(f (r), S(r)j) dr. 
Jt iev t 

Under Assumptions [231 this system has a unique solution, which we denote as (S l (i + s; t, Z t , St), £(i + 
s;t,Zt,St) ; s € R+,i £ Vt). We will denote the coordinate of the flow corresponding to ([5]) for an 
individual with initial condition £ as S(i + .; i, £, Zt, <Sf), so that in particular, S l (t + s; t, Z t , St) = S(t + 
s; t, Zt, St). In the sequel, we will also use that, if no births or deaths occur, for all < s < t: 

Z(t;s, H(s;0, €o,Z ,S ), Z t ,S t ) = S(t; 0, Z , So). (6) 

After a birth or a death the process is restarted with appropriately adapted initial conditions at that 
instant. 

The above description suggests a simple direct way for simulating the IBM. Starting from a birth 
or death event, first generate a standard exponentially distributed random number r, and then run the 
differential equations for the states of all Daphnia individuals and for the algae. Simultaneously integrate 
the sum of the birth and death rates of the Daphnia, starting from zero. When this integral reaches 
r, one of the Daphnia dies or gives birth. Which individual is the culprit and whether the event is a 
birth or a death is then decided from a single multinomial draw with probabilities proportional to the 
contributions of all the different events to the total event rate at that time. 

For the SDE that describes the process (Z(t), S(t)) t ^M + , we proceed as in [30], following a construction 
introduced by [T3] for the case without age or size, while accounting for the additional difficulty that 
the growth rate now depends on the rest of the population. To this end we again use that between 
two birth or death events the evolution of the population, conditionally on its state at the last event, 
is deterministic. As the integral form in which we present the SDE looks back at the end result of the 
events happening over a time interval, we also have to look in a retrospective manner at the resetting of 
the initial conditions at the moments that a birth or death occurs. Assume that the initial condition Zq 
and the initial food concentration Sq are given. The idea is that to construct the population at time t, 
we can proceed as follows: 

• If no event happens during [0, t] , then it is sufficient to consider the predicted sizes (S 4 (t; 0, Zq, So) ; i S 
Vo) of individuals at t. The population at time t is 

Zt=2_j %E«(t;0,.Zo,S'o),a*+t) 

»evb 

and the food concentration is St = E(t; 0, Zq, So). 

• If a birth event occurs at time s G [0,t], then V s = V s _ U {/ s _ + 1}, where I s _ is the number 
of labels already used so that the new individual gets labelled with the first available number 
j = I s _ + 1. The predicted sizes at time t, (E l (t; s, Z s _, S s _) ; i £ V s _), that we had for the 
individuals i £ V s _ before the event on the supposition that it were not to occur, are replaced by 
(S l (£; s, Z s _ +#(j,£ ,o)> S a _) ; i € V s _); see Fig. [TJ Moreover for the new individual with label j, we 
add a Dirac mass at (t; s, Z s _ + <5(j,f ,o)> S s _), t — s). 

• If the individual j £ V s _ dies at time s £ [0,t], then V s = V s _ \ {j} and the predicted sizes at 
time t, (H l (t; s, Z s _ , S s _) ; i £ V s ), that we had for the individuals i £ V s before the event on the 
supposition that it were not to occur, are replaced by s, Z s _ — 5q ^ a i yS s _) ', i £ V s ); see 

Fig. EJ Additionally, we delete the Dirac mass at (j, S J (t; s, Z s _ , S s _), a 3 s + (t~s)) that corresponds 
to the predicted size and age at time t of the dead individual. 
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Figure 1: At time there is a single particle of size £' and age a! , which is expected to have size £1 = 
H(£; 0, Zo, So) at time t, where Zo = 5(i,£' i0 ')- At time si, a second particle is born. Just before birth, the 
population is Z S1 _ = S^i.s( 31 :o,e ,z .s a ),a'+s 1 )- After the birth, the size expected for the first particle at time t is 
changed from £i to £2 = si, E(si; 0, Zo, So), Z S1 _ + 5^,( ,o), S S1 -) since there will be less resources for the 
two particles. At time S2, the second particle dies. The size expected for the first particle at t is changed again 
from £2 to £3 as there is now more resources for the first particles 's growth. 



Definition 2.2. Let Q(ds,di,d8) be a Poisson point measure (PPM) on M + x N* x K + with intensity 
ds Cg) n(di) <S> dO where ds and d9 are Lebesgue measures and where n(di) is the counting measure on N*. 
The PPM provides possible times of events. For each time, we draw the label i of the individual who 
may reproduce or die. The parameter 9 allows to define whether a birth of death occurs. Assume also 
that the initial condition are Z , So and Vo = {1,. . . , Io}, with E((Z n , 1)) < +00 a.s. Then 

^ = S (i,B'(t;0,Z ,S ),a' +t) + [ [ Q(ds,di,d0) % eVs 

(V s _+l,S Js - +1 (t;^Z s _+5 (/s +li5o , 0) ,S s ),t- S ) 

+ ( S {j^Ht;s,Zs_+S (Is _ +1 , (o , ),S e ),ai + (t-s)) ~ S (j,~.J(t:s,Z s _,S s ),ai + (t-s)))) 



where: 



x 1, 



fl<mi(v_,Z,_ ,S a ) 



+ ( - 6 (i,Z*(t;s,Z s _,S s ),ai + (t-s)) 

( S (j,E.i(t;s,Z s _-S {iti i_ ta i_ r S s ),ai+(t-s)) ~ (t;s,Z s _ ,S s ),aj+(t-s)))) 



(7) 



jev e _ 



x 1L 



,i(i,s_,Z,_ ,S s )<e<m 2 (i,s-,Z a _ ,S 3 ) 



mi (i,S-,Z s _,S s )=p(C_,ai_,S s _) 

m 2 (i,S-,Z s _ ) S s ) =m 1 (i,S-,Z s _,S s _) + fi(£ s ,a l s _, S s _). 



(8) 
(9) 



In this definition the first of the two terms in the square brackets corresponds to the births and the second 
to the deaths. The first term starts with the addition of one new individual, followed by the corresponding 
updating of the future course of the trajectories of the other individuals in order to eventually get the 
right outcome at the final time t. The first part of the second term, dealing with the deaths, starts with 
removing an individual, followed by a corresponding updating of the growth trajectories 
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Moment estimates obtained from ([7]) are very important for the proofs. By adapting the proofs in 
[2] (Th. 2.5) (see also [H]) where the main ingredient is the boundedness of the birth rate, we can show 
that: 

Lemma 2.3. Under A ssunwtions \2.1\ and i/E((Zo, 1)) < +oo as in Def. \2.2l then: 

(i) For anyT> 0, E( sup te[0>T] <Z t , 1)) < E«Z , l))e /3T < +oo, 

(ii) If we have additionally that E((Zo, l) p ) < +oo, then 

E( sup (Z f ,l) p ) < +oo. (10) 
te[o,T] 

As a consequence we obtain that (see [21 129] ): 

Proposition 2.4. The process (Z t , <S'()tgR + is well and uniquely defined on [0, T] for any initial condition 
such that E((Zq, 1)) < +oo and for any Poisson point measure Q. Moreover, the infinitesimal generator 
o/(Z t ) teR+ is flU). 



2.3 Martingale problem 

We end this section with a martingale problem that will be useful to derive the large population limits. 
Heuristically, the decomposition of the process (Z t ,f), for any test function / e C 1 , into a predictable 
finite variation process and a square integrable martingale can be viewed as a description of the paths 
as solutions of the evolution equation associated with the generator L (predictable finite variation part) 
plus noise (martingale part). The proof is given in Appendix IA1 

Proposition 2.5. Let us assume that E((Zo, l) p ) < +oo for p > 2. Let us consider a test function 
/(£,£, a) of class C . Then: 

M( ={Z t ,f(t, ., .)) - (Z 0) /(0, ., .)) - jf* (g( S) £ a) + yL(a, £, a) + 3 (£, £, a ) 

+ /(*, 6>, 0)/3(£, a, S.) - /(s, £, a)/z(£, a, S.)) Z.(d£, (11) 
is a square integrable martingale with predictable quadratic variation: 

(M f )t=J Q jf, (f 2 (s,Zo,0)f3(Z,a,S s ) + f 2 (s,Z,a)^,a,S s ))Z s (d(,da)d S . (12) 



3 Large populations 

We now focus on large populations of Daphnia. To that end we consider a sequence of processes describing 
the evolution of the population through time when starting with an initial condition of size proportional 
to the integer parameter K Q that we let increase to infinity. As already indicated, we may think of 
this K as the volume in which the population and its food live. We moreover scale the population with 
i.e., we transform from population size to population density. We thus consider a sequence Z K of 
populations such that: 



Z?(dj,d£,da) = — V S (l t hal) (dj,d£,da) (13) 



iGV/ 



where V t K is the set of individuals alive at time t. Again, Zf{di, da) = Zf (N* xd£x da) is the marginal 
on of Zj . We also consider (S^)t£R + the sequence of associated food concentrations, also indexed 
by K and satisfying: 

f(S?)~ i(Z,SK)Z t K (dj,dt,da), (14) 



dSl 

dt 

- ., „„ + 

with initial conditions (Sq)k&n* that converge in probability to Sq S M+. 

For each ifeN", KZ K has the dynamics of the process introduced in Dcfinition l2.2l with initial conditions 
KZq for which we assume that 



sup 

KeN 



^E((Z,f,l) 2 ) <+co. (15) 



G 



Proposition 3.1. Let T > 0. Under the Assumptions^^ the sequence [Z , S ) n >i introduced in this 
section converges in probability in D([0, T], A4p(M.^_) x ) to the unique continuous solution (C, g) of the 
following deterministic equations, characterized for any function f(t,£,, a) of class C 1 by: 

It tu \\ ir fin W-l. fir 3 ^ s "-) i df(s,.,.) d/(a,.,.) \ 
(Cuf(t, , .)> = (Co, f(0, , .)> + / {Cs, + + 9^^)ds 



+ / (Cs,f(s, ft, 0)/3(., .,g s )- f(s, ., ., e fl )>ds. (16) 
Jo 

§=/(*?(*))- / 7(e,e(*))Ct(de,do). (17) 



eft 

The proof of Proposition 13.11 is given in Appendix [Bj 

We conclude by showing that equations (fl~6]) - (fl~7]) allow us to recover the equations of [10]. We first 
establish a precise form of the solution £ t . That this measure is not absolutely continuous with respect 
to the Lebesgue measure on R^_, even if the initial condition Co is, was already noticed in e.g. ^Z§\ 130]. 
There it was stated that age and size where both parameterized by time, but no precise form for the 
measure for these cases was given. 

Corollary 3.2. Assume that the initial condition has a marginal measure in age that is absolutely con- 
tinuous with respect to the Lebesgue measure on K + so that C.o(d£,da) — qo(a,d^)da, where qo(a,d^) 
is a transition measure on K + . We denote by v a ,tid^) the image measure of qo(a — t,d£) through the 
application £ h-> H(t; 0, £, Co, £?o)- F° r an V t £ the marginal Ct at time t of the solution of p^|) is a.s. 
given by: 

Ct{d£,da) =la<t b(t - a)T{a,t - a, g[o,t])Ss{t;t-a,o,it- a ,et-a)( d O da 

+ i a >t F(t,a-t,£\Q[ 0)t ])lz =S ( t . £i iCo>eo )V ai t(d(3 da ( 18 ) 

where b(t) — L 2 f3(£,a, Qt)Ct(d£,da) is the total birth rate at time t, where Q[o t t] — (fo)se[0,*l an d where 
for a < t, 

F(a,t , Q[ 0l to+a]) = ex P ( ^ J Ai( s (io + a;t ,^ ,Ct , Qt ),a, g to+a )daj (19) 

is the probability that an individual born at to survives until age a when the food environment is given by 
(Qs)se[o,t +a]- For a>t, 



F'(t, oq, Q[ 0it] ) = exp ( - J ju(S(s; 0, f, Co, Qa), a + s, g s )ds^j 



(20) 



is the probability that an individual alive at t — with age ao and size £' survives until time t in an 
environment Q\o,t] ■ 

With Corollary 13.21 we recover the equations of [TUJ. Equation ^ provides the deterministic differ- 
ential equation describing the growth of Daphnias, represented by the distribution ( t (d£,da) (see ([18])): 

^(a) = g(Z(a), g (to + a)); £(0) = Co- 

If we consider an individual born at time to > and follow the survival probability through time a \— > 
T(a,to, £>[o,t +a]), Equation (JT9j gives the decay of the survival probability of an individual of age a at 
time t: ^ 

— (a, to, f?[o,t +a]) = — m(^(*o + «; to, Co, Ct , St ), a , £'to+a)-7 r (a, *o, £>[o,t +a])- 
From (|18p. we have the Daphnia population birth rate at time t: 

b(t) = I (3(E(t;t-a,0,Ct-a,Qt^a),a,gt)b(t-a)J 7 (a,t-a,g[ ,t])da 

/3(H(t; 0, C, Co, Qo), a + t, o t ).F'(t, a, £, 0[o,t])Co(d£> da) 
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The first term represents the contributions of individuals born after time 0, while the second term 
corresponds to individuals who where present initially. From this, we can deduce the algal concentration: 



dg t* 

= f(g{t)) - J 7( s (*; t-a,0, Ct-a, Qt-a), g(t))b{t - a)T{a, t - a, Q[ ,t])da 

7( s (i; 0, £, Co, Qo), Q{t))^'(t, a , (, Q[o,t])Co(d£, da). 



Proof, of Corollary \3. 6 A For the proof, we start by showing that £t(d4,da) admits a density m(£, a, i) 
w.r.t. a dominating measure underlying (|18l) . The equations satisfied by m(£, a, t) are then derived by 
separating the domain into R + x {a > t} and R+ x {a < t}, which corresponds to first studying the 
individuals born before and after initial time. 

First, recall that there is a unique solution g to (fT7|> . 
Let ip G C\ (R^_ , R) , let t 6 R+ and consider the associated test function: 

f(s,£,a) = (p(E(t;s,£,( s ,g s ),a + t- s). 

This function / is the unique solution of: 

+ («, £, «) = 0, f(t, e, a) = a) 

(e.g. 12 J. As a consequence, using this test function / in (fT5|): 

(C*,v) = (Co,^(S(t;0,e,Co,eo),. + *)> 

+ I" (v 5 ( s ( t ;s I Co,C s ,£ , s ),i-s)(G,/3(-,-,£' s )) 

s, £ i Csj &>), a + t - s)fj,(£, a, Q s )( s (d£, dafj ds. (21) 

The first term is related to individuals that are alive at time 0. The second integral relates to births 
between time and time t. The third term corresponds to the deaths between time and time t. 
If we consider positive functions <p>, then, neglecting the non-positive terms in (|21[) : 



0<(Ci,^)</ (/ tfi(Z(t; Q,£,(o, go), a)q Q {a-t,dO) da 

+ / b(t - a)ip(E(t;t - a,£, ,(t-a, gt-a),a)da. (22) 



Notice that the population is naturally divided into two sets. Since the aging velocity is 1, the individuals 
who were alive at initial time are of age greater than t at time t. Individuals born after time are of age 
smaller than t. So, if the function tp has support included in the set R+ x {a < t}, then we see that on 
R+ x {a < t}, Ct(d£,da) is absolutely continuous with respect to <fe(t;t-a,£o,ft-o,et-<0 (^4) da. Similarly, 
on the set R + x {a > t}, ( t (d£,,da) admits a density with respect to v a ,t{d£) da. Denote by m(£, a,t) the 
density of Ct with respect to the measure l a <t Ss(t;t-a4o,Ct-a,et-a)(^0 da + ^a>tv a ,t(dt;) da. 
Substituting this density in the third term of (|21[) gives, for ip with support in R+ x {a > t}: 

(Ct,<p) = / <p(E(t;0,£,(o,Qo),a)q (a-t,d£)da 



+ 00 

da I q (a-t, d£) tp(E(t;0,£, Co, go), a) 

+ 
t 

A*(S(s; 0, £, Co, Qo), a - t + s), g s ) 
m(S(s; 0, £, Co, Qo), a-t + s,s) 



ds 



By identification, the density m(£,,a,t) of Ct, on {a > t}, w.r.t. v a ^{d^)da satisfies: 



m(£,a,t)=l- / ii(E(s;0,£',( 0) Q ),a-t + s),g s ) 



n 



m(S(s; 0, Co, go), a - 1 + s, s) ds % =S (t ; o,e,Co,eo)} 



o 



where there exists a unique £' such that £ = S(i; 0, Coj £>o) under Assumptions 12. ll Notice that a — t 
is the age of the individual at time and is not a real function of time. Thus, we recognize an ordinary 
differential equation of degree 1 for s i— > to(S(s; 0, Coj go), a — t + s, s) from which 

m(£,a,t) =exp ( ~ ^ m(S(s;0,£',Co, go),a-t + s, s )ds) %=s(t;0,£',Co,eo)}- 

This yields the second part of (|T5|) . 

Choosing 99 with support in M + x {a < £}, (|21[) yields: 

(Ct,v)= / da b(t - a)ip(E(t;t - a,^o,Ct^a, et-a),a) 

- I da <p(E(t;t - a,£o,(t-a,Pt-a),a) 

fj,{E(u + 1 - a; t - a, £ , (t-a, Pt-a),u, p u+t - a ) 
m(S(w + £ - a; < - a, £ , Ct-a, Pt-a), w, w + t - a) 
da b(t - a)ip(E(t; t - a, Co, Ct-a, £>t-a), a) 

-/ da <p(E(t;t - a,£ ,(t-a, pt-a),a) 
Jo 

x / / fJ>(€,1J>, Pu+t-a) (u+t-a(d£,du) 
J JR+ 

Thus, on {a < t}, the density m(C,a,i) of Ct w.r.t. S^ t . t _ a ^ Q t _ a et _ a j(d^) da satisfies 
m(£,a,t)=b(t-a)- / u, p«+t-a)C«+t-a(rfC', dw). 



x 







du 











Notice that t — a is the time of birth of the individual and is not a real function of age. Then we recognize 
again an ordinary differential equation of order 1 from which 

m(£ , a, t) =b{t - a) exp (~ J p.(E(u + t - a;t - a, £ , Ct-a, Pt-a), p u +t-a)du^j . 

This ends the proof of the announced result (fT8|) for a < t. □ 

4 Conclusion 

In this paper we rigorously underpin the long suspected individual-based nature of physiological struc- 
tured population models as studied by the team of researchers surrounding Odo Diekmann. Such results 
were already known for purely age-based models (e.g. (27J [29l [30l H3 El HB] ) , and our present result still 
applies only to a subclass of the models considered by Odo and his co-workers, although a paradigmatic 
one, to wit their so-called "Daphnia" models e.g. [TUHl]. The i-state variables of these models move in a 
continuous deterministic fashion, dependent on their own value and the environment (food: "algae"), and 
births occur to a single fixed birth state at a rate that depends on the i-state and the environment (and 
are not e.g. coupled to specific i-state transitions). Even for this restricted class of physiologically struc- 
tured population models some twiddling of the existing probabilistic toolbox was in order. Yet, overall 
the toolbox proved its mettle, and the long standing assumption of a solid individual-based foundation 
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of the theory of physiologically structured populations was duly vindicated. The longer term goal is to 
extend this vindication to the full class of models put forward in [MJ El [5] . 

Acknowledgements: This work benefitted from the support from the "Chaire Modelisation Mathematique 
et Biodiversite of Veolia Environnement - Ecole Polytechnique - Museum National d'Histoire Naturelle - 
Fondation X" . 



Appendices 



A Proof of Proposition 12.5 



Let f(t,£,a) be a function of class C 1 . From ([7]), we obtain 
(Z t J) = f(t,Z i (t;0,Z< h S o ),ai+t) 

i&Vo 



Q(ds, di, d6) fljgy, 
i_n,y"-- .z .s.) (/(*> s (*5 s > °» + <*(/._ +u„,o), s s ), t - s) 

-f(t,Ei(t;s,Z s _,S s ),ai + (t-s))) 

+ ^ mi (i,s- ,Z B _ ,S s )<6<m 2 (i,s- ,Z S _ ,S B )[ ~ /(*> S ! %s- j Ss)> K + (* _ s )) 



■ (f(t,Z j (t;s,Z s _-6 (ut i0 j ): S s ),al + (t-s)) 
iev s _ 



f(t,^(t;s,Z s _,S s ),ai + (t-s)) 



(23) 



Using we have for any s < t: 



f(t, 3*(t; s, Z s , S s ), a* + (t - a)) =f(a, Q,<) + J (|{ + S'(«5 *, ^ <4 + « - a) 

+g(^ l (u; s, Z S ,S S ), S u ) — (u,^(u; s, Z s , S s ),a l s + u- s)Jdu 

Recall that we denoted by T k , k > 1 the birth and death events in the population. By convention, we let 
To = 0. Let us consider an individual i. Let to £ {Tc, k > 0} be the birth time of the individual (or if 
the individual is alive at time 0) and a\ be its age at time to (0 if to is the birth time). The sum of the 
terms in the r.h.s. of (1231) associated with individual i is equal to: 



_ /-tAT fc + 1 Vt c,r or 

/(*>.&.<)+E / (^ + ^ (s ' sl(s;Tfc '^' 5Tj 

k>0 JtAT k vt Oa 



)»«*„ +s-to) 



df \ 
■ g(s,^(s;T kl Z Tk ,ST k ),S s ) — (s,^(s;T k7 Z Tk ,S Tk ),a l to +s~t )jds 



l}'=i;i£V„_ ^mi{i,s-,Z B _ ,S 3 )<e<m 2 (i,s-,Z e _ ,S„)f( s i £s_ i a D <^2- 
JM*xK + J 
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The last integral correspond to the death term when individual i is dead before t. Thus, (|2"B1 gives: 



iev 



■ £ + _(., fl T k , Z Tk , S Tk ), a + s) 

k>0 tAlk 



■ - df ■ - \ 

■ g{s,^ l {s;T k , Z Tk ,S Tk ), S s)q^ ( s : - {s;T k , Z Tk , S Tk ), a l Q + s)j ds 

Q(ds, dj, d6)l j=i l l£V3 _ ]l mi(4jS _ : z 3 _ ,s s )<e<,n 2 (i,s-,z s _ ,s 3 )/0> £_ > a L ) 



JN*xl + 



(/(*,eo,o)+53 / 



JN*xI + 
t 

Q(ds,di,d0) l ieVs ^\ Vo 

-„-, (^ + f a ^- + \%T k ,Z Tk ,S Tk ),u-s) 
fc>0 tAJfcVs 

7 1 9 f I 1 \ 

+g{u, 2 /s - +1 (u; T k , Z Tk , S Tfc ) , S u ) ^ (u, S /s - +1 (u; T fe , Z Tfe , S Tk ), u - s) J du 

%< mi (i,s-,Z s _,S s ) 
~ /( S ) £s_ ) fl )1mi(!,s_ ,Z S _ ,S s )<e<m 2 (i,s_,Z s _ ,S S ) 

where the first bracket corresponds to individuals alive at time and where the second bracket correspond 
to individuals born after time 0. For s < u 



S (S'(u;s,z B ,s s )Mi+u-s)(d^da) = Z u (d£,da) 



iev s _ 

provided there has been no jumps between s and u. Thus, we have 

V .Qf Qf 



(Z t , fit, ., .)) =(Z , /(0, .,.))+ I 'da I Z s (d£, da) + £, a) + gfa S s M(s, £, a)) 



(/(*, Co, 0)/3(C, a, 5.) - f(s, £, a)ii(t, a, S.)) Z s (d£, da) 



1 



iev s 



f(8,£ ,Q)h<mi 



JN*xl + 

/( s , £s_ ) °)lmi(!,s_,Z s _ ,S,)<8<m2(i,s-,Z s _ ,S. 



(i,s-,Z s _ ,S B ) 

Q(ds, di, d9), 



where Q(ds, di, dO) = Q(ds, di, d6) — ds ®n(di) <g> d9 is the compensated Poisson point measure associated 
with Q. The integral with respect to (w.r.t.) Q(ds,di,d9) provides the martingale . This achieves 
the proof. □ 



B Sketch of the proof of Proposition 13.1 



When starting from (jllj) and using controls of moments as in [14] , the proof is similar to the one in [301 [22] • 
Step 1 We start by noticing that under the Assumption (|T5j) . we have the following estimate (e.g. [5]): 



sup fi( sup (Zl l ,l)'- 

nSN* Me[0,T] 



< +oo. 



(24) 



Moreover, from Assumptions 12 . 1 1 and ([S]), the size of any individual is bounded on [0,T] by £ — £o + gT 
and there exists for every e a non random constant S £ such that: 



sup p( sup S? > S s ) < e. 
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(25) 



From these estimates and Assumption 12 . 1 1 (ii) . there exists a constant A e € (0,A) such that: 

sup p( sup Z t n ([£o,£] x [0,A e ]) >e)<e. (26) 

nSN* V tS[0,T] ' 

Step 2 It is easy to see that the limiting values of (Z n , S ,n ) ne p}« are necessary continuous. Let us check 
the C-tightness (e.g. [IS]) of (Z n , S n )„ e ®* in D([0, T},M F (M. 2 + ) x M+). Using a criterion by 221 and given 
the compact containment that follows from Step 1, it is sufficient to prove the tightness of (S n ) ne w* and 
of the predictable finite variation part and martingale part of ((Z n , /))„ e pj* for / in C^(M^_,K) (which 
contains the constant function equal to 1). This is obtained by using Aldous-Rebolledo criteria (e.g. [19] ) 
and adapting for instance 2, 30 with the estimates of Step 1. 

Step 3 The identification of the martingale problem satisfied by the limiting values provides (Tl6T) - (fT7|) . 
Uniqueness of the solution of (|TB| - pT)) stems from the Assumptions 12.11 As a consequence, there is 
a unique limiting value and we have convergence in distribution of (Z n , S' n ) ra6 N* to the solution (£,q). 
Since the latter is deterministic, the convergence is also a convergence in distribution. 
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